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Abstract. - In this letter, we report progress on the field theory of polymerized tethered membranes. For 
the toy-model of a manifold repelled by a single point, we are able to sum the perturbation expansion in 
the strength go of the interaction exactly in the limit of internal dimension D — » 2. This exact solution is 
the starting point for an expansion in 2 — D, which aims at connecting to the well studied case of polymers 
(D — 1). We here give results to order (2 — D) 4 , where again all orders in go are resummed. This is a 
first step towards a more complete solution of the self-avoiding manifold problem, which might also prove 
valuable for polymers. 



Introduction. - The statistical mechanics of fluctuating lines and surfaces is a subject of great in- 
terest, which poses fundamental problems and has remained challenging for more than 20 years. One 
particular universality class, which has been studied extensively in the past, are polymerized or "teth- 
ered" membranes These are two-dimensional networks, where the bond-length 
fluctuates, but never breaks up. In the high-temperature regime nearest-neighbor interactions can be 
modeled by a harmonic potential. Neglecting self-avoidance, the membrane is extremely crumpled 
and highly folded, a property, which is characterized by the universal radius-of-gyration exponent v, 
defined as 

Rg ~ L" , V = , (1) 

where R g denotes the radius of gyration, and L is the linear internal size. Physically, < v < 1, but 
in the absence of interactions, the radius of gyration grows only logarithmically with the internal size. 

For a more realistic description one has to take into account self-avoidance, whose continuum ver- 
sion can be modeled by the generalized Edwards-Hamiltonian [|l0|] with 2-particle contact interaction 

H[r] = \ J (Vr(x)) 2 + | J J 8 d (r( X )-r(y)), (2) 

x£M xeMyGM 

where x£A4cM. D labels points in the manifold M, while r(x)GR d points to their position in external 
space. The Edwards model successfully describes long polymers [fQLO]. Much effort has been spent 
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to extend these results to membranes (D = 2). The problem is, that the usual e-expansion about the 
upper critical dimension is not feasible, since the latter is infinity. An important idea was therefore 
to generalize (Q) to manifolds of arbitrary internal dimension D. One then studies the D-dimensional 
manifold problem, and finally continuous analytically to D = 2. A major breakthrough was the 
proof of perturbative renormalizability |^|,^] to all orders in perturbation theory. This procedure was 
carried out to two loops [ |l3| , |l4] ] resulting in a radius-of-gyration exponent of v rj 0.86. This is a strong 
correction over the non-interacting theory with v = 0, but may still be in contradiction to Monte-Carlo 
simulations, which often but not consistently find tethered membranes in a flat phase with v = 1 JT^ , 
[T(|[T7|,[l8|]. While simulations are very demanding and therefore not yet conclusive, it is nevertheless 
compelling to try to identify possible mechanisms, which might render flexible membranes flat at all 
scales. Such a mechanism has indeed been found for rigid membranes, where fluctuations strongly 
renormalize rigidity Jl9|,[i[]. 

Here we study a simplified model, and solve it exactly at D = 2. It corresponds to a gaussian 
elastic manifold interacting by excluded volume with a single <5-like impurity in external space 



H[r\ = \ J (Vr(x)) 2 + g J S d (r(x)) . 



(3) 

xEM x£M 

As a first step to prove renormalizability of the full problem, [|| analysed (|j) and indeed showed 
renormalizability to all orders in perturbation theory for all dimensions < D < 2. (Q) has essen- 
tial features in common with SAM: Its critical embedding dimension tends to infinity as the internal 
dimension approaches D = 2. This can be read off from the dimension of the coupling go, which is 

[go] =:e = D- • (4) 

Thus, calculating universal quantities within the e-expansion necessitates similar techniques as for 
SAM, and we expect to learn more from the solution of the toy-model (||). 

Recently, we have been able | ]2l| ] to sum the perturbative expansion exactly in D — 2. The key-idea 
was, that when approaching D = 2, the correlator which enters all perturbative calculations, becomes 
essentially flat. In order to check the consistency of the results obtained by that method, one would 
like to go away from D = 2, and hopefully smoothly connect to polymers in D = 1, which are well 
enough studied to check almost any quantity. In [^T|], we have done a first step in that direction, and 
obtained quite promising results in first order in (2 — D), However, the expansion in (2 — D) is not 
a loop-expansion, and at each order in (2 — D), we have to resum an infinite number of diagrams. 
It turns out, that the results thus become very sensitive to the regularization procedure. In this letter, 
we pursue this road further, calculating contributions to the partition-function exactly for a manifold 
of toroidal or spherical shape. We obtain the expansion up to order (2 — D) 4 . This information can 
then be used to extrapolate away from D = 2. However, since we find that at D = 2, the fixed point 
is at infinity, one needs additional constraints, i.e. a scaling function, in order to be able to use this 
result. We have not been able to settle this question, despite the tremendous information contained in 
the perturbative result. We thus present our "raw data", together with some possible scaling-functions, 
encouraging the reader to think himself about the missing link. 

Perturbation theory. - Physical observables are derived from the partition function Z( go). We 
use it to define the effective coupling of the problem, 

g(z) := ^-(Z(0) - Z(g )) , (5) 
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which only depends on the dimensionless combination z := goL £ . Vm denotes the total internal 
volume of the manifold. Accordingly, the perturbation expansion reads 

re 00 / \N IN+l „ \ 

^^E^(ni/v(,,,) o , 

where the normalization of the (5-distribution has been chosen to be 5 d (r(x)) — (4Tr) d 8(r(x)) 
= J k c lkr ( x ') with J k := n~ d / 2 J d d k. Performing the averages within the gaussian theory with nor- 
malization y!_ /s d (r(x))) = 1, one arrives at 




N + l 



where C(x) :— ^ ({r(x) — r(0)) 2 ) denotes the correlator, and the S d ( J ^2 i ki) stems from the inte- 



gration over the global translation. Performing the shift k^+i — * fcw+i — 5Zi=i afl d integrating out 
the momenta k\, . . . , k^+i one obtains 

00 

l( z ) = Z Y] 777-TTj ( TT / ](detS))~ d / 2 , (8) 




(N+l) 

where we have factored out L e from the loop integration (such that the integrals now run over a torus 
of size 1), and the matrix elements are = ^[C(xN + i—Xi)+C(x]^ +1 ~Xj)~C(x i —Xj)]. 

Complete resummation of the perturbation series in D = 2. — Let us compute the iV-loop order 
of (Q): The asymptotic behavior of the propagator C{x) for large arguments is of the form 

C(x)~c + ±- In-, (9) 

Z7T a 

where cq denotes some positive constant (note C(x) > 0), and the logarithmic growth (for large x) 
is universal. In D = 2 we need an additional short distance cutoff a. The loop integrals, denoted 
by In, only depend on the dimensionless combination L/a. We can (somehow arbitrary) decompose 
det £ = (n£Li det S> with 



C(x N+1 -Xj) - C(xi-Xj) 



> 2 ' ^ 



C(xjv+i-Xi) 

1 . (10) 


(detS))- d/2 =: I N (L/a) = I? {L/a) (detl) (0) r d/2 . (11) 



The matrix 2> (0) denotes the limit a -» of @. It can be written as 3< ) = ±(I + ATP), where I 
denotes the identity and P the projector onto (1,1,..., 1), whose image has dimension 1, such that 
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det£)(°) = kpir-. Furthermore, to one loop I\{L/a) a =° ci(ln— ) d / 2 , where c\ denotes some 
(finite) constant. One then arrives at 

9{Z)=Z )^ Nl(l+N)*/™ ■ (12) 

A factor c\2 d / 2 has been absorbed into a rescaling of both z and g. The above series can be analysed 
in the strong coupling limit z — > oo. For this purpose we define functions f£(z) together with then- 
integral representation 



hW - 2 2_, ^(fc+Tv)^ ~ r(|) ' 



7V=0 



Thus in the limit of large z, the effective coupling ( |l2[ ) approaches the asymptotic form 



(14) 



Observables. - It immediately follows from this behavior that the correction-to-scaling exponent 
oj, which is defined as the slope of the RG-/3-function at the fixed point, equals zero. Here, it is useful 
to study the /3-function as a function of the bare coupling z, which reads (3(z) = — e z dg(z)/dz. 
Then, the correction-to-scaling exponent is obtained from the limit z — > oo of 

ez d[3{z) 
p(z) oz 

The value of u> can be checked in a Monte-Carlo experiment by considering plaquette-density func- 
tions on a membrane with self-avoidance in only a single <5-like defect. Be the partition function 
Z° = J V[r]S d (r(y)) exp[—H[r]], then the plaquette-density at the defect is obtained from (n) = 
dgfdz 'Sz^dz )' wnere ff = . O ne furthermore needs the density-density correlation at this point, 
which is defined as (n 2 ) ^ = gg/g z gi^Cff )■ ^ n me limit of strong coupling = ^(1+7) an d 
(n 2 ) = ^-(2+—+^), such that the ratio 



n )o «-°°. / e+u „,=o. ./I (16) 



V 2 £ +^ V 2 

becomes universal and should be measurable in simulations. 

(2 — D)-expansion. - Let us now analyse the theory below D = 2. Due to the renormalizability 
in < D < 2 and the existence of an e-expansion we expect the renormalized coupling to reach a 
finite fixed point in the strong coupling limit as soon as D < 2. This approach is characterized by a 
powerlaw decay of the form 

g (z) = g* + S(\nz) + Oiz-^/ 6 ) , (17) 

where S is some scaling-function growing at most sub-exponentially and u>\ > uj > 0, with co defined 
in (|l5|). In order to gain information about g below D = 2 one has to expand the loop integrand 
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(det D) d l 2 in powers of 2—D. For convenience, we take a — > 0. The propagator takes in infinite 
D-space the form C(x) = \x\ 2 ~ D / (S D {2-D)), where S D = 2ir D / 2 /T{§) denotes the volume of the 
Z?-dimensional unit-sphere. The factor (Sd(2— D))^ 1 replaces ln(-j) and is absorbed into a rescaling 
of the field and the coupling according to r — > r (So (2—D)) and go - * So (So (2— D)) d / 2 , such 
that the factors of (In ^)~ d/2 in (Q and @ are replaced by (S D (2-D)) d / 2 . The propagator in the 
rescaled variable can then be written as 



C(x) = 1 + (2 - D) C(x) 



(18) 



where for convenience of notation we allow C(x) to depend itself on D. 

Of course, on a closed manifold of finite size, C(x) needs to be modified, but the form ( |l8[ ) is 
independent of the shape of the manifold. Accordingly, one may expand the matrix D, which is 
53 = + (2—D) U>, where is defined as before and coincides with the limit D^2 when 
inserting the above C(x) into 2). Moreover, D is of the same form as D, but each C(x) has been 
replaced with C(x): = b[C(xjv+i — ^i) + C(xat + i — Xj) + C(xj — a^-)]. Then, 



detS = detl) (0) exp|Tr ln(l + (2 - D)^' !]-^) | 



(19) 



where [S)' )] -1 = 2(1— -^-j-P) denotes the inverse matrix of f)(°\ Expanding the integrand ( |l9| ) in 
powers of (2 — D) and the coupling g , all orders in g can again be summed, with the difference 
that the integrands are no longer constant. Expanding up to the nth order in 2—D involves n powers 



of C(x). Introducing the notation f(xi, . 



,Xk) 



Xl, 



,Xk) with the integration 



defined as f := J d D x (on the torus) the overbar can be thought of as an averaging procedure. To 
first and second order in 2—D, the only integrals to be evaluated are C(x) and C 2 (x). In order to 
reveal the structure of the expansion we generated all terms up to fourth order. Generally, the terms 
are of the following form 



]T(det£> 



(Oh-d/2 



N=l 



nLi (Tr([l)(0)]-iB)".)"Vz) 
(N+Tjl 



E 



ft +2j (z) 



m m 2 ■ • • raj • 
il »2 ■ ■ ■ "I , 



ff +2j (z), (20) 



"1 m 2 ■ ■ ■ m; max 
l l "2 ■ ■ ■ "i | min 



where max and min are some integers, and summation over the index j is implicit. The precise form 
of the vector entries MP will be reported elsewhere [22]. The renormalized coupling then reads up to 
fourth order in 2—D (note that we have absorbed a factor of 2 d l 2 in both g and z): 



g(z) = /^(z)-(2-D)-M^ |;) /r 2J 0) 
+ (2-D) 2 



-{2-Df 
-(2-Df 



0(2 - L>) 5 (21) 
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From the integral representation ( 13 ) of ff +J (z) and the above expansion, it follows immediately that 
the exact renormalized coupling can be written as 



g(z) = z dr g{r) e~ 



where g(r) is of the form 

9{r) = 



T{ d±2 ) 



n—Q j— — n mELX 



(22) 



(23) 



Let us try to gain more information about the powerlaw behavior in (|17|), that is about the expansion 
in 2— D of the correction-to-scaling exponent lo. Powerlaw behavior forces the series (^3|) to turn into 
some exponentially decaying function g(r) as can be seen from the asymptotic form of g(z) 



oo 



r(i+f) 



(24) 



Now, we test a possible form of the exact g(r), which is consistent with the expansion (|Tj) and which 
satisfies the following properties: (i) In the limit of D — 2 the exact form r d / 2 /r(^t^) emerges and 
(ii) for D < 2 the corresponding g(z) has a finite fixed-point value together with a strong coupling 
expansion. The (non-unique) ansatz is 



;(r)=C 



1 -£(£>, r) e~ 

uo/e 



d/2 



(25) 



where S(D, r) is analytic in D = 2 of the form S{D, r) = 1 + |r 5„(r)(2-D) n , and each 

S n (r) has a Laurent expansion S n (r) = s «/ 



Note, that in the limit of D 



2, the 



expression (g5|) gives r d / 2 , while for D < 2 it yields upon integration the form (|4|), ensuring both 
properties (i) and (ii). Inserting uj/e = LU2{2—D) 2 +0(2—D) 3 (the linear term in (2— D) has to 
vanish) into the ansatz (g5b and expanding to second order in 2— D provides 



<J 



(r) = C r d / 2 



1-- [S l {r)[2-D)+ 



U>2 



d-2 



S 1 (r) 2 +S 2 (r))(2-D) 2 + 



(26) 



The first coefficients of the (2— D)-expansion of g(r) obtained from (eTJ) read 



Kr) 



r d/2 



r( d±2) 



l + (2-£))|c (!-£-) - (2 D) 



2r 



d 2 (d 



2^ r 



C 2 +xC r~ 2 



(27) 



Comparing @ and @, one identifies C = l/r(^), 5i = -C(l-|i) and to 2 = 2C 2 , where 

C c (x):=C(a;) — C. Note that the terms proportional to C in iSaC?*) mostly cancel with 5i(r) 2 , a sign 
that the ansatz catches some structure. 

The diagrams to be calculated at this order are C and C 2 . On a manifold of toroidal shape, which 
is equivalent to periodic boundary conditions, two discrete sums have to be evaluated: 

1 



= Sr 



y- 



fc#0 



2tt(2-L>) 



-0.44956 + 0.3583 (2 - D) + 0{2 - D) 1 



C 2 



S dJ2^4 = 0.152661 + 0(2 - D) . 



(28) 
(29) 
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k is /J-dimensional with components hi = 2tt/L rii, and m integer. With the results given above, this 
leads to 

to = 2eCf (2 - Df + 0(2 - Df = 0.305322 e (2 - D) 2 + 0(2 - Df , (30) 

which can be compared to the exact result for D = 1 (polymers): co = e. As a caveat, note that the 
above scheme is not unambiguous, since different ansatze in (^5|) are possible. Also the second order 
term proportional to r in ( p7| ) could in principle either be attributed to lu 2 or iS 2 - More constraints are 
necessary to settle this question. 

In summary: We have presented a complementary approach to treat the problem of tethered mem- 
branes in interaction. We hope that this approach will prove fruitful for self-avoiding tethered mem- 
branes, with eventual applications for polymers. 

* * * 
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